Local C 0,a Estimates for Viscosity Solutions 
of Neumann-type Boundary Value Problems 
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Abstract 

In this article, we prove the local C 0,a regularity and provide C°' a estimates for 
viscosity solutions of fully nonlinear, possibly degenerate, elliptic equations asso- 
ciated to linear or nonlinear Neumann type boundary conditions. The interest of 
these results comes from the fact that they are indeed regularity results (and not 
only a priori estimates), from the generality of the equations and boundary condi- 
tions we are able to handle and the possible degeneracy of the equations we are able 
to take in account in the case of linear boundary conditions. 
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1 Introduction 

In this article, we are interested in the local C 0,a regularity of viscosity solutions of 
nonlinear Neumann boundary value problems of the form 

F(x,u, Du, D 2 u) = in O, , , 

G(x, u, Du) = on 90, ^ 
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where O C R n is a smooth domain, F and G are, at least, real- valued continuous functions 
defined respectively on O x R x R n x S n and 90 xlx R n , where 5 n denotes the space of 
real, n x n, symmetric matrices. The solution u is a scalar function and Dm, D 2 u denote 
respectively its gradient and Hessian matrix. More precise assumptions on 0,F,G are 
given later on. 

We recall that the boundary condition G = is said to be a nonlinear Neumann 
boundary condition if the function G satisfies the following conditions 

(Gl) For all R > 0, there exists /i R > such that, for every (x, u,p) G dO x [-R, R] x R n , 
and A > 0, we have 

G(x, u,p + Xn(x)) — G(x, u,p) > /i R X , (2) 

where n(x) denotes the unit outward normal vector to dO at x G dO. 

(G2) For all R > there is a constant K R > such that, for all x,y G (90, p,q£ M n , 
u,v E [—R,R], we have 

|G(x,u,p)-G(y,u,g)| <iirfl[(l + b| + |g|)|a;-?/| + b-g| + |u-u|] . (3) 

The main examples of boundary conditions we have in mind are the following : first, 
linear type boundary conditions like oblique derivative boundary conditions, in which G 
is given by 

G(x,u,p) = {p,-f(x))+g(x), (4) 
where 7 : dO — > W n is a bounded, Lipschitz continuous vector field such that 

(~/(x),n(x)) > f3 > for all x G 90, 

and 5 G O 0,/3 ((9O) for some < f3 < 1. Here and below, "(p, g)" denotes the usual scalar 
product of the vectors p and g of W l . 

Next nonlinear boundary conditions : the first example is capillarity type boundary 
conditions for which G is given by 

G(x,u,p) = {p,n(x))-e(x)y/T+W, (5) 

where 8: dO — > W 1 is a Lipschitz scalar function, such that \9{x)\ < 1 for every x G dO. 
A second example is the boundary condition arising in the optimal control of processes 
with reflection when there is control on the reflection, namely 

G(x,u,p) = sup{(>y a (x),p) + c a (x)u- g a (x)}, (6) 

aeA 

where A is a compact metric space, 7 a : dO — > M™ are Lipschitz continuous vector fields 
such that ( i y a (x),n(x)) > (3 > for all x G dO, and c a ,g a : dO — > R are Lipschitz 
continuous, scalar functions. 
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We are going to show that, under suitable assumptions on O, F, G, any continuous 
viscosity solution of §B) is in C^{0) for some < a < 1, with an estimate on the 
local C°' Q -norm of u. These results are indeed regularity results (and not only a priori 
estimates); this is their first main advantage. But they are also valid, in the case of 
linear boundary conditions, for possibly degenerate equations, a second original feature. 
The counterpart is that the regularity properties we have to impose on O, F, and G, are 
stronger than in the case of a priori estimates where the solution u is assumed to be either 
in C 2 (0) n C^O) or at least in W 2 ' n (0) n C l (p). 

The classical a priori estimates for this type of problems are in fact proved for linear 
equations and extended to fully nonlinear equations by a simple linearization procedure 
described in Lieberman & Trudinger [16]; this linearization requires the regularity of the 
solution. For linear equations, the classical results are obtained under rather weak assump- 
tions on the coefficients of the operators in the equation and in the boundary condition. 
To the best of our knowledge, the first results in this direction are the ones of Nadirashvili 
[201 12T] for linear, uniformly elliptic equations with L°° coefficients associated to oblique 
derivative boundary condition of the form 

(Du,<y(x)} + a(x)u(x) + g(x) = in 30. (7) 

He first proves them with a continuous direction of reflection 7 and for Lipschitz domains, 
and then for a direction of reflection in L°° for C 2 domains. Results in this direction were 
also obtained in the 80's by Liberman [H] by different methods. Recent improvements 
on the regularity of the coefficients of the equation (which can be assumed to be L p for p 
large enough or in L n ), were obtained by Kenig & Nadirashvili [IB] and Lieberman [15j . 

The case of fully nonlinear equations was first considered in Lions & Trudinger [19] who 
show the existence of a smooth classical solution in C 2 (0) nC 1,1 (0) for Hamilton- Jacob i- 
Bellman equations with smooth coefficients and directions of reflection. As mentioned 
above, in Lieberman & Trudinger [16], the case of fully nonlinear equations is considered 
in a more systematic way but most of the results are obtained by a linearization procedure 
and are based on results for linear equations; it is worth pointing out that, in general, 
the passage from a priori estimates to regularity results requires the existence of smooth 
enough solutions for a sequence of approximate problems (and even uniqueness for the 
problem itself), and is often only valid for convex or concave equations. 

Our approach is based on classical viscosity solutions (not L p viscosity solutions): for 
a detailed presentation of the theory of viscosity solutions and of the boundary conditions 
in the viscosity sense, we refer the reader to the "Users'guide" of Crandall, Ishii and Lions 
[8] and the book of Fleming and Soner [9], while the books of Bardi and Capuzzo Dolcetta 
[T] and Barles [2] provide an introduction to the theory in the case of first-order equations. 

Clearly, this approach requires more regularity properties for the operator G which 
has to be locally Lipschitz with respect to its variables, for the domain O which has to 
be assumed to be C 2 and for the equation (F has to be continuous). Its advantage 
is that, in the case of oblique boundary conditions of the form (J7|), we require only 



3 



F(x,u,p, M) to be nondegenerate (in a sense precised below) in one direction which 
depends, near the boundary, only on p and 7. Whereas in the case of more general 
nonlinear boundary conditions we require F to be uniformly elliptic. To prove such 
results, we use systematically an idea introduced in Ishii & Lions [12] which has already 
been used to obtain interior, local regularity (or global regularity) in [3] and Barles & 
Souganidis [6]. 

In this paper, for technical reasons, we treat separately the "linear case" , i.e. typically 
the case of oblique derivative boundary condition where the operator G is linear with 
respect to u and p and the "nonlinear" case where G is not linear. A surprising fact in the 
linear case - and maybe our result is not optimal in this direction - is that the assumptions 
on F, and in particular the ellipticity one, depends on 7. We were unable to remove this 
dependence. 

In the case of Neumann boundary condition i.e. when j(x) = n(x), our "strong 
ellipticity condition" can be written formally as 



where A > and where, here and below, the notation p stands for — . This condition 

\P\ 

is the natural requirement for the interior C 0,a regularity to hold and it allows to extend 
the results up to the boundary. 

Next, if the direction of reflection 7 is C 2 , then a classical property which is used in 
Lions [T7] (see also Lions & Sznitman [18J) is the existence of a C 2 function A(x), taking 
values in the set of non-negative symmetric matrices and such that A(x) r y(x) = n(x) for 
any x G dO. In this case, we have to require that the above "strong ellipticity condition" 
is valid but replacing in (EJ) p by A^ 1 (x)p. Since A is not unique, this assumption is 
admittedly not completely satisfactory. 

Finally, if j(x) 7^ n(x) is just Lipschitz continuous, then we have again to assume 

"strong ellipticity condition" of F in the A~ x {x)p direction where again A(x) r y(x) = n(x) 
but, here, A(x) is just Lipschitz continuous and this creates technical difficulties.. 

We mention that, both in the linear and nonlinear case, we prove the regularity result 
by assuming that G does not depend on u. Indeed one can always reduce to this case by 
a suitable change of variable that we show later on. 

The proofs of these results rely on the constructions of suitable test-functions inspired 
by the test-functions built for proving uniqueness results : in the case of Neumann or 
regular oblique derivatives problems, the corresponding uniqueness results were proved by 
Lions [17] (see also [8]) and in the case of nonlinear Neumann-type boundary boundary 
condition in [3]. It is worth pointing out anyway that the construction in the case of 
Lipschitz continuous 7's, which is the difficult case, takes a completely different form 
here. 

It is worth mentioning also the results of Ishii [11] proved, in the case of nonlinear 
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Neumann boundary conditions, under weaker assumptions on O but stronger assumptions 
on the boundary condition than in [1]; our approach requires more regularity of the 
boundary and therefore we do no use the test-function built in [11] . 

This paper is organized as follows : in Section 2, we state our regularity results both 
in the linear (Subsection 2.1) and in the nonlinear case (Subsection 2.2) and we provide 
the main proofs. Such proofs rely on the constructions of a suitable test-functions which 
are different in the linear and nonlinear case: these constructions are given in Section [3j 
It is worth pointing out anyway that, despite most of the arguments are common in these 
two cases, the conclusion is a little bit different because of the particular " ellipticity 
conditions" used in these two cases. 

2 The Local C° " Estimates 

In this section, we state and prove the local C°' a regularity of the solutions of the problem 
(JT]) both in the linear and nonlinear cases. As pointed out in the Introduction, these two 
cases requires slightly different assumptions. We first introduce the assumptions which 
are common of both cases. First, for the domain O, we require 

(HI) (Regularity of the boundary) O is a domain with a C 2 -boundary. 

This assumption on O implies the existence of an IP-neighborhood V of dO such that 
the signed distance function d which is positive in O and negative in O c is in C 2 (V). We 
still denote by d a C 2 -extension of the signed distance function to R n which agrees with 
d in V and we use below the notation n(x) = —Dd(x) even if x is not on the boundary. 

The "strong ellipticity" conditions on F are different in the linear and nonlinear cases 
but the following natural growth condition on F is, on the contrary, the same 

(H2) (Growth Condition on F) For any R > 0, there exist positive constants Cf, 
Cg, Cf and functions , ug, zu R :R + —>R such that ujf (0+) = 0, uj§{r) = 0(r) 
as r ^ +oo, zu R (t) — > as t — > +oo, and for any x, y E O, —R < u,v < R, 
p, q e W 1 , M e S n and K > 

F(x,u,p,M) - F(y,v,q,M + Kid) < cjf (\x - y\(l + \p\ + \q\) 

+ ™ R (\P\ AM)|p-g|)||M|| 

+UJ *(K)+C? + C?(\p\ 2 + \q\ 2 ) 
+C*\x - y\(\p\ 3 + \q\ 3 ) , 

where \p\ A |g| = min(|p|, |g|). 

In the sequel, K always denotes a positive constant which may vary from line to line, 
depends only on the data of the problem and is, in particular, independent on the small 
parameters we are going to introduce. 
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2.1 The case of linear boundary conditions 

In this subsection we examine the case when G is linear with respect to p, namely it is of 
the form ([7j). 

The main additional assumptions on F and G are the following. 

(H3a) Oblique-derivative boundary condition and ellipticity : there exists a 
Lipschitz continuous function A: O —>■ S n with A > Cold, for some cq > such 
that A(x)j(x) = n(x) for every x G dO, and for any R > 0, there exist L R , X R 
such that, for all iGO, \u\ < R, \p\ > L R and M,N G S n with iV > 0, we have 

F(x,u,p,M + iV) -F(x,u,p,M) < -\ R (NA^~(x)p, A^~(x)p) + o(l)||JV|| , (9) 

where o(l) denotes a function of the real variable \p\ which converges to as \p\ 
tends to infinity. 

Before stating the assumption on the boundary condition, we want to point out that 
the existence of such A is really an assumption in a neighborhood of dO, then, under 
suitable assumptions on F, A can be extended to O and even to M. n . 

For the boundary condition , we require 

(H4) (Regularity of the boundary condition) The functions 7 and a in (J7D are 

Lipschitz continuous on dO, ( i y(x),n(x)) > [3 > for any x G dO and g is is in 
C%£(dO) for some < j3 < 1. 

Our result is the 

Theorem 2.1 Assume (Hl)-(H2)-(H3a)-(H4) . Then every continuous viscosity solu- 
tion u of (TJP with G given by (0) is in Cf^(O) for any 0<a<lif(3 = l and with 
a = j3 if (3 < 1. Moreover the -norms of u depend only on O, F, 7, a, g through 
the constants and functions appearing in (H2)-(H3a), the local C 0,1 -norm of 7 and a, 
the local C '^ -norm of g and the local C 2 -norm of the distance function of the boundary 
including the modulus of continuity of D 2 d. 

Despite we make a point here to have a rather unified result as we do it for the proof, 
Theorem 12.11 contains clearly three cases which are rather different from the technical 
point of view. 

1. The homogeneous Neumann boundary condition — — = or more generally — — \-g(x) = 

on on 
when g is a C 2 -function, is the simplest case. Of course, one can take A as being 

the identity matrix and the construction of the test-function does not require the heavy 

regularization procedures we use in Section El 

It is worth pointing out here that, in the construction of the test-function, the two 
du 

terms — and g are treated in fact separately. To prove a result with g being either Holder 
on 
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or Lipschitz continuous requires the rather sophisticated regularization argument of case 
(3) and is therefore of a different level of difficulty. 

2. The case of "regular" oblique derivative boundary condition does not differ so much 
technically from the first case. The assumption which says that A(x)j(x) = n(x) where, 
for any x G dO, A(x) is a nonnegative symmetric matrix with a C 2 dependence in x, 
implies that 7 is a C 1 function of x. 

3. The case when 7 is only Lipschitz continuous and when g is only Holder or Lipschitz 
continuous, is technically very different as it is for the comparison results (cf. [1]). Here 
the only (known) way to treat this case is through a non-trivial regularization argument 
which we are going to use also. 

Of course, in the proof below, we emphasize more the (difficult) third case : the proofs 
are far easier in the two first ones. 

Assumption (H2) is a classical hypothesis in such type of regularity result : it is the 
same as the one which appears for the interior regularity (cf. Ishii & Lions [12], Barles[3j). 
It is worth pointing out, anyway, that the treatment of the oblique derivative boundary 
condition does not lead to a stronger assumption. 

Concerning (H3a), we recall that, for the interior regularity, just the strong ellipticity 
"in the gradient direction" (cf. (jHJ)) is needed as in the case of homogeneous boundary 
condition. Unfortunately, in the case of oblique derivatives boundary conditions, this 
natural assumption does not seem to be enough or, at least, it has to be reformulated 
in a far less natural way. Of course, all these conditions hold when a classical uniform 
ellipticity property holds, like (H3b) below. 

One of the main examples we have in mind is the case of standard quasilinear equations 

-Tt[b(x,Du)D 2 u) + H(x,u,Du) = in O, (10) 

where b is a n x n matrix and H a continuous function. In this case, the assumptions 
(H2) and (H3a) are easily checkable. 

(H3a) is equivalent to : there exists A > such that, for any x G O, p G M n 

b(x,p) > \A~ l (x)p <g> A~ 1 (x)p - o(l) Id , 

where, as in (H3a), o(l) is a function of \p\ which converges to as \p\ — ► +00. 

This assumption may seem restrictive, in particular the fact that the constant A does 
not depend on x and p; but, for general equations, if b satisfies the above statement with 
a strictly positive A depending on x and p, one can divide the equation (i.e b and H) by 
X(x,p) and the above property becomes true with A = 1. 

We conclude these remarks about (H3a) by emphasizing the role of the "o(l)" term 
and, for the sake of simplicity, we assume that A = Id. Without this term, (H3a) would 
be essentially reduced to 

b{x,p) > \p® p , 
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for any x G and p G IR n — {0}, while, with this term, (H3a) is satisfied if 

b(x,p) > Xq(x,p) (g) q{x,p) , 

where g is a continuous function such that \p\~ 1 {q{x,p) — p) — > as p — > oo, uniformly 
with respect to x G O. This condition is not only a more general assumption on b but it 
is also far easier to check it. Of course, a similar remark can be made for general A's. 
Now we turn to (H2). It is satisfied when 

(i) b is a bounded, continuous function of x and p and there exists a modulus of 
continuity uj\ : IR + — > M + and a function w : 1R + — > R + such that w(t) — > as £ — > +oo 
and 

|6(x,p) - < Wiflz — + |p| + \q\) + w(\p\ A |g|)|p- q\) . 

Moreover the uniform bound on b provides a c^2 with a linear growth. 

(ii) The function H satisfies : for any R > 0, there exist positive constants C± , Of, 
Cjf such that, for any x,y G O, —R < u,v < R and p, q G M n , 

f/>, it, p) - H(y, v, q) < C* + C?(\p\ 2 + \q\ 2 ) + C*\x - y\ (|p| 3 + |g| 3 ) . 

As we already mentioned it above, this assumption is classical (See Ishii & Lions [T2]. 
Barles[3]). 

Proof of Theorem 12.11 : We are going to do the proof in two steps : the first one 
consists in proving the result when a = and contains the main arguments. In the second 
one, we show how to handle the u a(x)u" term. According to the proof of the comparison 
result, it is clear that to take in account such a w-term in regularity result is not immediate 
and we are going to do it in a very indirect way, by adding an extra variable. 

Step 1 : The a = case. 

Since we are going to argue locally, we start with some notations. For every x G O 
and r > 0, Bq(xo,t) := B(xo,r) fl O denotes the open ball in the topology of O while 
Bq(x , r) := B(xq, r) fl O denotes the closed ball in the topology of O and 8Bq(xq, r) : = 
dB(x , r) n O. In a similar way we define Bq{x , r), B o (x , r), and dB o (x , r). 

We assume that g G C^(dO) for some < /3 < 1, the case g G Cf^(dO) being treated 
in a similar way. 

We are going to prove that, if we choose a = j3, for all x G O, if r is small enough, 
then there exists a constant C depending on the different data of the problem, such that, 
for any x G Bq(x , r), we have 

u(xq) — u(x) < C\x — x \ a . (11) 

All together, these inequalities give the answer provided that we control the dependence 
of r and C in xq, which will be the case. 

The proof of this estimate is done in two steps whose arguments are the same: the 
first step consists in proving that the result holds for a small enough (depending on the 
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local L°° norm of u and the data of the problem) and then that this property implies that 
the result holds for a = (3. 

We provide the main arguments of the proof of (fTTj) in the case when the boundary 
condition plays a role, i.e. when B(xo,r) D dO ^ 0; the other case is simpler and can be 
treated by the methods of In the common proof of the two steps, we therefore argue 
with some a < (3. The following lemma is the key stone of the proof. 

Lemma 2.1 Assume that B(xq, r) (IdO ^ and that u is a bounded, continuous solution 
of (QJj in Bq^xq^t) with the oblique derivative boundary condition on B(x ,3r) n dO. 
Under the assumptions of Theorem \2.1\ on F , 7 and g, then there exists a constant C > 0, 
depending on F,j,g and \ \u\\L°°(B—( Xo ,3r)) suc ^ that for all x G Bq(xq,t) the estimate 07]) 
holds. 

Proof. In order to prove (II ip . we consider the auxiliary function 

$o(x, y) = u(x) - u{y) - o (x, y) , 
where the function Bo has the following form 

e (x,y) = Ce-^ +d ^[M^yP + Le- R ^ +d ^M^xo)+Xo(x,y) , 

where a G {0,(3} is a fixed constant, C, L, K are some large constants to be chosen later 
on and where the continuous functions ip (x,y), Xo( x ,y) satisfy the properties listed in 
the lemma below. In order to point out the main dependences in these functions but also 
to simplify the rather technical estimates we have to make in the proofs, we introduce the 
following notations which are used in all the sequel 

X := X ^ - , Y := x — y , Z :— d(x) — d(y) , T := d(x) + d(y) . 

Lemma 2.2 Under the assumptions of Theorem \2.1\ for 5 > small enough, there exist 
real-valued, continuous functions ip$(X,Y,T) , x$(X,Y,T, Z) defined respectively in M n x 
E n x R and m I" x I" x I x R, such that if we set 

ip s {x, y) ■= ^s{—^—, x-y, d{x) + d(y)) , 

Xs(x, y) := Xs{-^— , x-y, d(x) + d(y), d(x) - d(y)) , 

the following facts hold for some constant K depending only on the local C 0,1 -norm 0/7 
and the C ^ -norm of g. 

(i) For any 1,761", T G R and for 5 > small enough, 

K-^Yl 2 < ip s (X,Y,T) < K\Y\ 2 + K5 , (12) 



9 



-K\Z\ < xs(X,Y,T,Z) < K\Z\ + K5 a . (13) 

(it) When 5 - 0, 4>s(X,Y,T) - $ (X,Y,T) and X s(X,Y,T,Z) -> x (X,Y,T,Z), 
uniformly on each compact subset o/K" x 1" x 1 and M n x 1" x 1 x 1 respectively. 

(Hi) When S > 0, the functions ips, Xs are C 2 functions and the following estimates 
hold for any X, Y, Z, T 

(D^MX, Y, T)Y, Y) = 2$ 5 + OQY | 3 ) + 0(5) as Y — , S -> 
(D Y i>s(X, Y, T),Y) = 24> s + 0(\Y\ 3 ) + 0(6) as Y — , 6 -> 

|D| X ^(X,F,T)| <K|Y| , |££ y ^(X,Y,T)| <tf|Y| 
|D yy ^(X,F,T)| < X, , \D T MX,Y,T)\ < K5 
Dtt^s = D^ x ip 5 = Dryips = , 

Moreover if A s := (\Y\ 2 + 5 2 )^ 

iDxUX^Y^Z^KA^lZl , 

\D z xs(X,Y,T,Z)\ < K , \D T xs(X,Y,T,Z)\ < K6, 
\D 2 X8(X,YT,Z)\<A^ 2 \Z\ . 

(iv) There exists a constant K > large enough (independent of C and L) such that, if 
we set 

Q s (x,y) = Ce-*W^»[V*(z,^ , (14) 

then, for \x — y\ small enough, we have 

(D x Q s (x,y), 1 (x))+g(x)0 if x e dO , (15) 
(-D y e s (x,y)^(y)) + g(y)<0 ifyedO. (16) 



The proof of the key Lemma [2.21 is postponed to Subsection 3.1. The first two prop- 
erties in the point (iii) in Lemma 12.21 are going to play a central role in the proof. 
We continue with the proof of Lemma 12.11 

We are going to show that, for a suitable choice of L > 0, chosen large enough in order 
to localize, then for C > large enough, we have 

M LyC :=_ max $ (x,y)<0 (17) 

Indeed if (IT7|) holds then plugging x = x and using the estimates (fl2"j) and ([TBI we get 

u(x ) - u(y) < K\x - y\ a , 
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for some constant K > depending on a, r, \\u\ \L^(S- u { Xo ,r))i F through (H2)-(H3a), 
the uniform Lipschitz norm of 7 and the uniform Holder norm of g in Bq(xq, r) and 
independent of x (at least if x remains in a compact subset of O), which is the desired 
regularity result. 

To prove (|T7|) . we first choose L, C > large enough in order to have $>o(x, y) < for 
x or y on dB(x ,r) R O. This is possible since u is locally bounded on O and since the 
conditions (|T2|) and (|T3|) imply that 

-2/| 2 < ^o(x j2 /) < C 2 |x-2/| 2 , (18) 

-C 3 \x-y\ < Xo(x,y) < C 4 \x-y\ . (19) 

Of course, L and C depends on r. 

From now on, we fix such an L and we argue by contradiction assuming that, for all 
C > 0, Ml^c > 0. Since $0 is a continuous function, the maximum is achieved at some 
(x, y) G B-Q(x ,r) x Bq(xq,t) and we observe that, by the choice of L, C, we may even 
assume that x G B-Q(x ,3r /4) and y G B-Q(x ,r). Here we have dropped the dependence 
of x, y on C for simplicity of notations. 

Two quantities are going to play a key role in the proof 

Q i: =C\x-y\ a , 
Q 2 := L\x - x \ 2 , 

(again we have dropped the dependence of Qi,Q 2 in C for the sake of simplicity of 
notations). The reason for that is the following: by using only the local boundedness of 
u, we are only able to show that Qi,Q 2 are uniformly bounded when C becomes very 
large while if we use the local modulus of continuity of u, we can show that Q%, Q 2 — > as 
C — > +00. The idea of the proof can therefore be described in the following way: we first 
show that u is locally in C 0,a for a small enough with suitable estimates depending only 
on the local L°° norm of u and on the data, and this is done by using only the uniform 
boundedness od Qi,Q 2 . Then this first step provides us with a local modulus of continuity 
for u and we obtain the full result using this time that Qx, Q 2 — > as C — > +00. 

As we just mention it, from the fact that $ (x,y) > 0, using classical arguments, 
Qi, Q 2 are bounded and, more precisely the following estimates hold, in which K denotes 

the constant f e 2Kr {2\ \u\ ^ + HxolU) J 

Q x < K a or equivalents \x-y\< KC~ 1/a , (20) 
L\x — x \ 2 < K . 

From the first estimates (1201) it follows, in particular, that \x — y\ —>■ as C — > +00. These 
estimates hold true for any maximum point of the function $0 m B-Q(xo,r) x Bq(xo,t). 
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The function Bo is not (a priori) a smooth function and therefore we cannot use 
directly viscosity solutions arguments; this is why we have to consider the functions ips 
and xs defined in Lemma 12.21 Since ipg — > ipo and x$ ~~ ¥ Xo as ^ ~~ *• locally uniformly in 
OxO, for all L, C > 0, there is 6~c,l > such that for < 5 < 5c,l, we have 

__ max (u(x) — u(y) — Qs(x, y)) > . (21) 

Let (x s , y s ) be the maximum point of the function (x, y) \— > u(x) — u(y) — Qs(x, y) 
in Bq{xq,t) x _Bq(xo,t). Standard arguments show that, up to subsequence, (x s ,y s ) 
converges to a maximum point (x, y) of $ as ^ "~ > 0. Moreover we may suppose that 
x s — y s 7^ 0. Indeed if for all 5 > we have x 5 — y 5 = 0, then x — y = as well. But in this 
case we would have <& Q (x,y) < which is a contradiction. Hence, we can assume without 
loss of generality that x s — y s remains bounded away from 0. For simplicity of notations 
we now drop the dependence of (x s , y s ) on 5 as we already dropped it on C. For C large 
enough, we have x, y £ Bq(xq, r). Moreover, from Lemma [2.21 it follows that 

(D x Q s (x, y),j(x)) + g(x) > if x e dO , (22) 
(-D y Q s (x, y), j(x)) + g(y) < if y G dO . (23) 

Thus the viscosity inequalities associated to the equation F = hold for u(x) and u(y) 
whenever x, y lie. 

By the arguments of User's Guide [8], for all £ > 0, there exist {p,B{) E J ' u(x), 
(<7> -S2) £ ' u(y) such that 

p = D x Q s (x, y) , g = -Z? J/ e 5 (a:, y) , 

-(e- 1 + || J D 2 e 5 (a;, 2/ )||)/d< ( ^ _^ ) < D 2 5 (a;, y) + e(D 2 Q s (x, y)f , (24) 

and 

F{x,u{x),p,B 1 )<0 , F{y,u{y),q,B 2 ) > 0. (25) 

We choose below e = p\\D 2 Qs(x,y)\\~ 1 for p small enough but fixed. Its size is 

determined in the proof below. Next we need the following lemma. We recall that Y 

Y 

denotes x — y and Y = rpr. 

Lemma 2.3 If p is small enough and if B\, B2 satisfy (Efy then, for \Y\ small enough 
(i.e. for C large enough), there is K > such that 

{{B l -B 2 )Y,Y) < -CK- 1 a(l-a)^- + 0(6)\Y\ a - 4 (26) 

+CK\Y\ a - 1 + K , (27) 
as 5 — * 0. Moreover Bi — B\ < Ki(Y, 5)1 d with Ki(Y, 5) given by 

K\iY, 5) = K{C5\Y\ a - 2 + C\Y\ a ~ l + 1) . 



12 



Proof of Lemma 12.31 By the regularity properties of xs given in Lemma 12.21 it 
is tedious (but easy) to check that all the terms in D 2 Qs(x,y) are estimated by K + 
Ci^|y| a_1 except perhaps the ones coming from the derivation of the first term. More 
precisely, we have 

D 2 Q s (x, y) = M 1 + M 2 + M 3 + M 4 

where 

M 1 = c e -*««H%)) [^sf^D^s + |(| - m s r^D^ s ® Dips 

M 2 = C^s) a/2 D\e- RWx)+d ^) , 
M 3 = Ca(ii 5 y /2 -\D(e- R ^ +d ^) ® Dips) , 
M 4 = LD 2 (e- Rd{x) M^ x )) + D 2 X6 (x, y) . 

But we have also 

D x ip s {x, y) = * 6 + Dy'ips + D T 4>sDd(x) , 

Dx ips ~ ~ 

D y ips(x, y) = — ^ h Dy ^ 5 + D T^sDd(y) , 

D 2 xx M*,y) = ^^ + D 2 YY ij 5 + D 2 XY ^ 5 + D T ^D 2 d(x), 

D 2 yy ip s (x,y) = ^^ + D 2 YY i> s -D 2 XY i> 5 + D T ^sD 2 d(x), 

D 2 xy M^V) = ^f^-Dyy^g. 



and, taking in account the properties given in Lemma 12.21 it can be readily checked that 
||M 2 || < CK\Y\ a } ||M 3 || < CKadYl " 1 + 5\Y\ a - 2 ) and ||M 4 || < K . 
On the other hand, for all £, ( G M n , we have 

(m 1 (c, o, (e, 0) = ce- k ^ +d ^ - o, (e - o) + we + o, « - 0) 

+(p 3 (e +o,(e+o + we, o, (e, o>] (28) 

where 

P l = D 2 YY ^ 5 ) a/2 (X,Y,T) , P 2 = 2D 2 XY (4, 5 y/ 2 {X,Y,T) , P 3 = D 2 xx ^ 5 ) a/2 (X,Y,T) , 

and P 4 is the matrix involving all the terms D 2 XT {ip s ) a / 2 , D YT (-ijj s ) a / 2 , D^ T (ijj s ) a ^ 2 to- 
gether with the -Dd and _D 2 cZ derivatives. One can easily check that P 4 is estimated 
by 0(5)(tps) ( ' a ~ 4: ^ 2 ■ On the other hand, Lemma [2.21 implies 



\Pi\ = 0(tys) 



(a-2)/2\ 
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\P*\ = 0((^s) ia - 1)/2 ) , 
\Ps\ = 0{^s) {a - 1)/2 ) ■ 
Choosing £ = ( in ( |28i) . we first deduce that 

M x < C (0(^ s ) {a ~ 1)/2 ) + 0(5)(^) (a " 4)/2 ) N ) Id . 



We next choose £ = — ( = Y. According to the two first properties in the point (iii) of 
Lemma [2.21 and taking in account the above estimate on P4, we get 



(M l (Y,-Y),(Y,-Y)) < 4C c - 2 *W+*» aW ' s) " /2 



D> Y j t Y,Y) + £-l)m^ 
1 4>s 



2 \Y\ 2 
+C(P A (Y,-Y),(Y,-Y)) 
< CK~ l a{a - l)^r 2)/2 + K(^s) a/2 + 0(6)^t i)/2 . 
By combining the above estimates, we obtain 
(B 1 - B 2 Y, Y) < CK~ l a{a - \)\Y\ a ~' 2 + CO{8)\Y\ a ^ + CK\Y\ a ~ l + K . (29) 

And the final upper estimate on Bi — B 2 follows from the estimates on Mj for i = 
1,2,3,4. □ 

We continue with the proof of Lemma 12.11 : we estimate \p — q\, \p\, \q\ and 
||£>i||, H-B2II. For some K > 0, we have 

\p\, \q\ > CK~ l a\Y\ a ~ l - 0(5)\Y\ a - 3 - K , 

\p\, \q\ < CKalYl*- 1 + 0(5)\Y\ a - 1 + CK\Y\ a + K , 
|p - q\ < KCa\Y\ a + 0(\x - x \) + o s (l), as \Y\ -> 0, 5^0, 

\\B 1 \\,\\B 2 \\<K{1 + -^--) [Ca\Y\ a - 2 + 0(6)\Y\ a ~ 4 + l] . 

At this point, it is worth noticing that we are going to let first 5 tends to for fixed C 
and we recall that, since we assume that M LjC > 0, Y does not converge to when S 
tends to for fixed C. The first consequence of this fact is that the term 0(<5)|y| a_3 is 
playing no role in the lower estimate of |p|, \q\ since we can choose 5 as small as necessary 
and, by the above estimates, we have \p\, \q\ — > +00 as C — > +00 

We are going to use (H2)-(H3a) with R = WuWloo^b^^^ . We drop the dependence 
in R in the coefficients and modulus which appear in these assumptions. We subtract the 
two inequalities (125]) and write the difference in the following way 

F(x, u(x),p, B x ) - F(x, u(x),p, B 2 + Ki(Y, 5)Id) 

< F(y, u(y), q, B 2 ) - F(x, u(x),p, B 2 + Ki(Y, 6) Id) , (30) 
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and, using the fact that B\ — B 2 < Ki{Y, S)Id, we apply (H3a) to the left-hand side and 
(H2) to the right-hand side of (1301) . Recalling also that \q\ — > +00 as C — > +00, this 
yields 

XTi[(B 2 -Bi + Ki(Y, 6)Id)(jf\x)p ® A^^xjp)} + o(l)\\B 2 -B x + Ki(Y, 5)Id\\ 
< u x {\x - + \p\ + \q\) + w(\p\ A \q\)\p - q\)\\B 2 \\ + u 2 {K 1 {Y 1 5)) + C x 

+ C 2 (\p\ 2 + \q\ 2 ) + C 3 \x - y\(\p\ 3 + \q\ 3 ) (31) 

Now we use the following result which is a consequence of the construction of the 
test-function and whose proof is given at the end of Subsection 3.1. 

Lemma 2.4 We have 

A r Hxjp = Y + oy(l) + 0,5(1) as \Y\ -> 0,5 -> 0. 

We want to point out that, in the above lemma, \Y\ — > is in fact equivalent to C 
going to infinity. 

We come back to (13~TT) : by Lemma [2.41 and recalling also that \p\, \q\ — > +00 as C — > 

+CX3, we get 

Tr[(5 2 - 5i + ^i(y, ^/^(A^p ® jf^{x)p)} 

> ((B 2 - B X )Y, Y) + K^Y, 5)-\\B 2 - ^11^(1) + o,(l)) . 

Moreover by using the estimates on Si , B 2 we have 

[Ca|F| Q - 2 + 0(5)|F| Q - 4 + 1] . (32) 

As we already pointed out above, we are going to let first 5 tends to for fixed C and, 
since we assume that M LjC > 0, Y does not converge to when S tends to for fixed C. 
Hence in the estimates below, we are going to replace the terms which converge to as 
5 — ► by 0,5(1). On the other hand, as we already mention it above, C going to infinity is 
equivalent to Y going to and when C is going to infinity, p, q are also going to infinity; 
we can therefore incorporate the o(l)-term coming from (H3a) in the oy(l) term. 

Therefore, by combining fl32l) . the estimates on ||-B 2 ||, \p\, \q\ and \p — q\, and Lemma 
12.41 we are lead to 

Tr[(5 2 - B x + K\(Y, 5)Id)(A^{x)p ® A^Uxjp)] 

> CK- l a{\ - a)\Y\ a - 2 + o 5 (l) - K 

- (CKa\Y\ a - 2 + 0,5(1) + K)(o Y (l) + o,(l)) . 



\B 2 — Sill < K 



1 + 



O(p) 
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On the other hand, for the right-hand side of (13 ip . we first look at the u\ term. By 
tedious but straightforward computations, we have 

\x — y|(l + |p| + \q\) + zu(\p\ A \q\)\p - q\ = KaQ 1 + Kzu(\p\ A \q\)Q^ 2 + o y (l) + o 5 (l) , 

1 /2 

since 0(|x — xq\) is like Q 2 • This estimate is emphasizing the role of Q\,Qi and the 
necessity of having the w term. 

The complete estimate of the right hand side of (!3~T|) is 

Kuji (KaQi + Kw{\p\ A \q\)Ql /2 + o Y (l) + o s (lj) Ca\Y\ a - 2 

+ KC 2 a 2 \Y\ 2a ~ 2 + C 3 a 3 |y| 3Q - 2 + K + o y (l) + o*(l) , 

where we (partially) use the fact that Q 1 = C\Y\ a is bounded for C large enough. 

By dividing all the above inequalities by the (very large) term Ca|y| a_2 , we obtain 
the following (almost) final estimate 

XK- 1 < Ku t (KaQi + Kw{\p\ A \q\)Q l J 2 + o y (l) + o*(l)) + KaQ 1 + Ka 2 Q\ 

+ o Y (l) + o s (l) . 

And by using the fact that \p\, \q\ — > +oo as C tends to +oo, this yields 

XK' 1 < Kuj x (KaQx + o Y {l)Q)l 2 + oy(l) + 05(1)) 

+KaQ 1 + Ka 2 Q 2 1 + o Y (l) + o s (l) ■ (33) 

Using this last estimate, the conclusions of the two steps we mention at the beginning 
of the proof follow rather easily. 

On one hand, by using the uniform control on Qi,Q2, we can choose a small enough 
(depending only on the local L°° norm of u and the data) in order to have 

^K- 1 > - (Ku^KaQi) + KaQi + Ka 2 Q\) > Kuj x {KaQ x )) + KaQi + Ka 2 Q\ . 

With this choice, it is clear that the above inequality cannot holds for 5 small and C large 
enough (depending again only on the local L°° norm of u and the data) and the local C 0,a 
estimate is proved for small enough a. 

On the other hand, repeating this proof for any x G l%(xo,2r), this C°' a property 
provides us with a modulus of continuity in Bq(xq,t) (which depends only on the L°° 
norm of u in Bq{xq, 3r) and the data), and in the above estimate, for any a < (3, we can 
use the fact that Qi, Q 2 — > as C — > +00. Arguing as above, we obtain the C°' a estimate 
for any a < (3. And the proof of the Step 1 is complete. □ 



16 



Step 2 : how to handle the ll a(x)u v term 

We are going to introduce an extra variable to reduce this case to the previous one. More 
precisely we consider the function d:Ox1-4K defined by 

v(x,y) = exp(ky)u(x) , 

where k > is a large constant to be chosen. 
This new function is a solution of 

3 2 v 

+ exp(ky)F(x, exp(— ky)v , exp(—ky)D x v, exp(—ky)Dl x v) + k 2 v = in O x R , 

and 

(7(2;), D x v) + -a(x)D y v + exp(ky)g(x) = on 80 x R . 

We first remark that since we are going to argue in a neighborhood of the point (x , 0), 
the exponential terms are not a problem to check the assumptions on either the equation 
or the boundary condition. The only difficulty comes from (H3a) since we have changed 
the boundary 7 and therefore A is not given anymore. 

In order to define the new matrix A, denoted below by A k since it depends on k, we 
first set n := (n, 0), the exterior unit normal vector to dO x R and 7 := (7,/c" 1 a). In 
fact, because of the form of (H3a), it is more convenient to define A^ 1 and we do it by 
setting 

r v / A~ 1 (x) fc" 1 a(a:)n(a:) \ 
Ak [X) — \ k~ l a(x)n T (x) 1 J ' 

where n T is the transpose of the column vector n. 

An easy computation shows that A^ 1 (a;)n(a;) = ^f(x). Moreover, if P = (p,p n+ i) E 
W l+1 , we have 

(A- 1 (x)P,P) = (A- 1 {x)p,p) +2k- 1 a{x)(n{x),p)p n+l + \p n+1 \ 2 , 

and applying Cauchy-Schwarz inequality to the second term of the right-hand side, it is 
straightforward to show that, for k large enough, A^ l {x) is still a definite positive matrix. 
Finally we consider the operator F:OxRxKx R n+1 x S n+1 — > R defined by 

F(x, y, v, P, M) = -M n+hn+1 + exp(ky)F(x, ex.p(-ky)v, exp(-ky)p, exp(-ky)M) + k 2 v , 

where M is the n x n symmetric matrix obtained from M by removing the last column 
and row. 

Next we claim that the operator F satisfies (H3a) with A k 1 but by replacing the 
op(l) term which tends to as \P\ — > 00 by a op(l) + 0^(1) where 0^(1) — > as k — > +00 . 

Indeed, we first observe that, since the new solution v is Lipschitz continuous w.r.t 
y and since y is a variable which corresponds to a tangent direction to the boundary, 
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the property ll \P\ — > +00" is equivalent in fact to u \p\ — > +00" because p n+ \ remains 
bounded. 

On an other hand, by using (H3a) for F, one can easily see that the checking of our 
property reduces to show that, for all R > 0, x G O, P = (p,p n+ i) and ( = (£, 77) G K n+1 , 
we have 

X R (A- l (x)p,0 2 + N 2 > \ R (A^(x)PX) 2 - o k (l) , 

for some constant A# > 0. Because of the particular form of A^ 1 (x), this property is 
obvious for u k = +00" and, of course, this implies that it is also satisfied for k large 
within a Ofc(l)-term. 

We finally observe that the norm of A^ 1 does not depend on k if we choose it 
large enough. In order to conclude, we just remark that the proof of the first step still 
works if the term o&(l) is small enough and the proof is complete. □ 



Remark 2.1 We remark that, under further regularity assumptions on O and the coef- 
ficients appearing in the boundary condition ([TJ), it is possible to handle the a(x)u term 
without adding an extra variable but by using another change of variable. More precisely, 
let us suppose that the following assumption holds 

(H5) There exists a C 2 -function x '■ O — > K such that — = a(x) on dO . 

cry 

Then the function v defined by v{x) = e x ^'- ), u(x), is a solution of a modified equation in 
O (but still satisfying (H2)-(H3a)) with the boundary condition 

— + e x{x) g(x) = on 90, 
(77 

and we can apply the proof of the Step 1 of Theorem 12. II to v. 

Assumption (H5) holds for example in the following case : if O is a C 2,/3 domain 
and 7, a are C 1,13 function, then the existence of x is given by Theorem 7.4 (p. 539) in 
Lieberman Sz Trudinger[16j. Indeed, to build x, one can solve the Laplace Equation in O 
together with the oblique derivative boundary condition. 

Moreover, if O is a C 3 domain and 7, a are C 2 , then one can just take 

. a(x)d(x) 
(n(x),j(x)) 

in O where a, 7 and n denotes here suitable extensions of these functions to O. 
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2.2 The case of nonlinear boundary conditions 

In this subsection, we consider the case of nonlinear boundary conditions of the form 

G(x,u,Du) = in 90, (34) 

where G: dO xlx R" — > R is a continuous function, satisfying the conditions (Gl) and 
(G2). 

In this section, we use the following assumptions on F and G. 

(H3b) (Uniform ellipticity) For any R > 0, there is X R > such that, for all x G O, 
-R <u < R, p G R n and M,N G <S n such that M < N, we have 

F(x, u, p, M) - F(x, u, p, N) > \ r Ty(N - M) . 

(G3) For all R > and M > there is K R}M > such that 

dG 

\(-Q-{x,u,p),p) - G(x,u,p)\ < K R>M , (35) 

for all x G dO and for all p G R n , \p\ > M, \u\ < R. 
(G4) There is a function G^ : 80 x R x R n -> R such that 

—G(x,u,\p)^G oc (x,u,p) as A — > oo . (36) 
A 

locally uniformly in (x, 

Before providing our result, we want to point out that the Goo appearing in (G4) is 
homogeneous of degree 1 and satisfies (Gl) and (G2). 
Our result is the 

Theorem 2.2 Assume (Hl)-(H2)-(H3b) and (G1)-(G4). Then every bounded con- 
tinuous solution u of (Hp is in C^"(0) for any < a < 1. Moreover the -norms of u 
depend only on O, F, G, through the constants and functions appearing in (H2)-(H3b) ; 
and in (G1)-(G4) ; the local C 2 -norm of the distance function of the boundary including 
the modulus of continuity of D 2 d. 

Proof of Theorem 12.21 We are going to do the proof in three steps : in the first one 
we prove the result in the case when G is independent of u and homogeneous of degree 1 
in p, then in the second step, we remove the homogeneity restriction and finally, in step 
3, we use the method of the second step of the proof of Theorem 12.11 to deal with the 
dependence in u. 

Step 1 : The case when G is independent of u and homogeneous of degree 1 
in p 
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The proof is similar to the one of Theorem 12.11 and we just outline the main differences. 
Again we treat only the case when the boundary plays a role. 

Since the boundary dO is C 2 , by making a suitable change of variables, we can assume 
without loss of generality that the boundary is flat and more precisely that OC\B(xq, 3r) C 
{x n 0} and dO fl B(xo,3r) C {x n = 0}. It is worth noticing that the assumptions made 
on F and G are preserved by such a change. In order to keep simple notations, we still 
denote by F, G the functions arising in the equation and in the boundary condition in 
the domain with flat boundary. 

We have to prove the following lemma which is the key stone of the proof. 

Lemma 2.5 Assume that B(xq,t) fl {x n > 0} ^ and that u is a continuous solution of 
(TJP in B(xq, 3r) fl {x n > 0} with nonlinear boundary condition G = on B(xq, 3r) fl {x n = 
0}. Under the assumptions of Theorem \2.2\ on F, G, there exists a constant C > 
depending onF,G, \ \u\\L°°(B—(x ,3r)) such that, for any x e B(xo, r)fl{x n > 0} ; the estimate 
(COP holds. 

Proof of Lemma 12.51 . In order to prove ffTTl) . we consider the auxiliary function 

&o(x, V) = u(x) - u(y) - o (x, y) , 
where the function 0o has in this case the following form 

e {x, y ) = Ce-^ x " + y^[Mx,yP + i4 {x,x ) , 

for some large constants C, L, K to be chosen later on and where the continuous functions 
ijjo(x,y), 4>o{x,y) satisfy the properties listed in the following lemma in which we use the 
notations 

v x + y v 

X := , Y := x — y . 

2 

We want also to point out that the parameters 5 and rj we introduce in this lemma 
play completely different roles in the proof, the role of S being far more important than 
the role of t] which is a small but fixed parameter; this is why we choose to drop the 
dependence in rj of the functions i/jg, 4>s below. 

Lemma 2.6 Under the assumptions of Theorem \2.2l there is a function <p e (^(IR 2 ™) 
and, for 5 > and r\ > small enough, there exists a real-valued, continuous function 
ips(X,Y) defined in R n x R n such that if we set 

ij) S (x, y) := x-y) 

the following facts hold 

(i) There exists a constant K > such that, for 5 > small enough, 

K-^Yl 2 < MX,Y) < K\Y\ 2 + K5 , (37) 
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for any Y = (Y 1 , . . . ,Y n ) el". 

(ii) When 5 — > 0, ips(X, Y) — > ^o(X, Y) uniformly on each compact subset ofW 1 x R n . 

fraj When 5 > 0, the function xps is C 2 and the following estimates are valid for some 
constant K > 

\DxMX,Y)\ <K\Y\ 2 + K5, \D Y MX,Y)\ < K\Y\ , 

(D Y MX, Y),Y) = 2MX, Y) + 0( V )0(\Y\ 2 ) + 0(6) , 
(D 2 YY MX, Y)Y, Y) = 2MX, Y) + 0( V )0(\Y\ 2 ) + 0(5) 
as Y — > 0, r] — > and 5 — > . 

\D 2 XX MX,Y)\ < -\Y\(6 2 + \Y\ 2 ) 1 ' 2 , |L>|^(X, Y)| < K\Y\ , |L>f^(X, Y)| < K 
V 

(iv) 

K- 1 \x-y\ A <M^y)<K\x-y\\ (38) 
ft^ There exists K > /arge enough (independent of C and L) such that, if we set 

Q s (x, y) = Ce-*^^ ty s {x, y)]§ + L<j> (x, x ) , (39) 

then, for \x — y\ small enough, we have 

G(x,D x e 5 (x,y))>0 ifx n = 0, (40) 
G(y,-D y e B (x,y))<Q ify n = Q. (41) 

The proof of the key Lemma [2761 is postponed to the Subsection 3.2. To prove (iTTj) . we 
first choose L, C > large enough in order to have $o(x, y) < for x or y on dB(x , r) DO. 
This is possible since because of the conditions (1371) and (138|) . Of course, L, C depends 
on r. 

As in the proof of Theorem 12.11 we fix such an L and we argue by contradiction 
assuming that, for all a G (0, 1) and C > 0, M^c > 0. Since $o is a continuous function, 
the maximum is achieved at some (x,y) G Bp(x ,r) x Bq(xq,t) and we observe that, by 
the choice of L,C, we may even assume that x G Bq(xo, 3r/4) and y G 5o(xo,r). Here 
we have dropped the dependence of x,y on C for simplicity of notations. 

We use here Q\ := C\x — y\ a and Q2 := L\x — xo\ 4 . From the fact that $o(x,y) > 0, 
using classical arguments, the following estimates follow, in which K is the constant 

/ \ l/a 

(2e 2 ^|| W ||oo) 

Qi < K Q or equivalents |x - y | < (42) 
Q 2 <K 
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These estimates hold true for any maximum point of the function $ i n Bq(xq, r) x 

Since the function $o is not (a priori) a smooth function we have to consider the 
functions defined in Lemma T2.6I The following property holds : for all L, C > 0, there 
is Sc,l > such that for < 5 < 5c, l and K large, we have 

max (u(x) — u(y) — Qs(x,y)) > . (43) 

Bg{iO,r)x%(iO,r) 

Let (x 5 , y 5 ) be the maximum point of u(x) — u(y) — Os(x, y) in Bq(xo, r) x Bq(xo, r). 
Standard arguments show that, up to subsequence, (x s , y s ) converges to a maximum point 
(x, y) of $o as 5 — > 0. Again we may suppose that x 5 — y s ^ 0. For simplicity of notations 
we drop the dependence of (x s ,y 5 ) on 5 as we already dropped it on C. 

For C large enough, we have x, y G B-Q(xo,r). Moreover, from Lemma [2761 it follows 
that 

G(x, D x Q s (x, y)) > if x n = 0, (44) 
G(y,-D y Q s (x,y))<0 i£y n = 0. (45) 

Thus the viscosity inequalities associated to the equation F = hold for u(x) and u(y) 
whenever x, y lie. 

By the arguments of User's Guide [8], for all e > 0, there exist (p, Pi) G J ' u(x), 

2 — 

(q, B 2 ) G J ' u(y) such that 

p = D x Q s (x, y) , g = -P^B^x, y) , 

-(e- 1 + || J D 2 6 5 (a;, 2/ )||)/d< ( ^ _° B2 ) ^^e^^+ep^^,!/)) 2 , (46) 

and 

F{x,u{x),P,B x ) <0 , F{y,u(y),q,B 2 )>0. (47) 

We choose below £ = p||D 2 9,5(x,?/)||" 1 for p small enough but fixed. Its size is 
determined in the proofs below. In order to have estimates on B\ and P2, we set $(X, Y) = 
e~ KXn (ips(X, Y))i, with the correspondence given above between X, Y and x,y. 

By the regularity properties of tps, <Po given in Lemma I2.6[ it is tedious (but easy) 
to check that all the terms in D 2 Qs{x,y) are bounded except perhaps the ones from 
D 2 ijj(X, Y). Therefore, the inequality (|4"6"|) can be rewritten as : for all £, ( G lR n 

(b^,o-(b 2 c,o < (i+o(p))[<Pi(£-c),(£-c)> + we+cMe-o) 

+(p^+o,(t+c))+m\ 2 + \<\ 2 )], (48) 



for some constant K and where 

Pi = L%y1>{X, Y), P 2 = D\ Y i>{X, Y) , P 3 = P| X ^(X, F) . 
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Moreover, as we remark above, we can assume that Y = x — y remains bounded away 
from and this implies (after again tedious computations) that P3 is bounded as well. 
Choosing £ = ( in the above inequality, we first deduce that 

B x -B 2 < K 2 (Y,5,rj)Id 

with 

K 2 (Y,5,ri) = K[C (a^Y^- 1 ^ 2 + |F| 2 ) 1/2 + S\Y\ a ~ 4 ) + l] . 

We next choose in (l4~8j) . £ = — ( = Y. By using the properties on the first and second 
derivatives of ips and O proved in Section 3.2, we get, for some K > 



((B 1 -B 2 )Y,Y) < KCe 



-2kx a{ifj s )' 1 



2 \Y\ 



1 Vs 



K 



< -KCa{\ - + CaO(r])\Y\ a ~ 2 + 0(5)\Y\ a - 6 + K. 

If (ej)i<j< n _i are {n — 1) vectors such that (ei, . . . , e n _i, Y) is an orthonormal basis of lR n , 
we know that 

n— 1 

Tr (B 1 - B 2 ) = J2(( B i ~ B ^ e i) + " ^2)^, f) , 

i=l 

and by combining the above estimates, we deduce 

Tr (Si-S 2 ) < -C , Ka(l-«)%^ + 0(5)|F| Q - 6 + C'aO(r ? )|r|^ 2 + ^ 2 (r,5,r / ) . (49) 

|r I 

Finally, for \Y\ small enough (i.e. for C large enough), we get 

Tr (B 1 -B 2 ) < -CKa(l - a)\Y\ a - 2 + 0{5)\Y\ a - 6 + CaO(7])\Y\ a - 2 + K 2 (Y,5, V ) . (50) 

Now by using the estimates on the first and second derivatives on 9,5 shown in the Sub- 
section 3.2, we get, for some K > 0, 

\p\,\q\ > Ca(K- 1 -Krj)\Y\ a - 1 + 0(6)\Y\ a - 3 -K, 
\p\,\q\ < CaK\Y\ a - 1 + 0(5)\Y\ a ~ 1 + K, 
\p-q\ < CaK\Y\ a + 0(5) + 0(\x - x \ 3 ) , 

||Si|U|S 2 || < K(l + -^)(l + Ca\Y\ a - 2 + 0(5)). 



As in the proof of Lemma 12. 1[ we notice that we are going to let first S tends to for 
fixed C and we recall that, since we assume that M L C > 0, Y does not converge to 
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when 5 tends to for fixed C. The first consequence of this fact is again that the term 
0(5)|y| a_3 is playing no role in the lower estimate of \p\, \q\ since we can choose 5 as small 
as necessary. The new point in the above estimate is the 7/-term : we choose it sufficiently 
small in order to have, say, K' 1 — Kr\ > . With this choice, by the above estimates, 
we have \p\, \q\ — > +00 as C — > +00 

We are going to use (H2)-(H3b) with R = 1 1 ii || 2,°° (B^o.r) • We drop the dependence 
in R in the coefficients and modulus which appear in these assumptions. We subtract the 
two inequalities (|4T|) and write the difference in the following way 

F(x, u(x),p, BJ - F(x, u(x),p, B 2 + K 2 (Y, 5, rj)Id) (51) 
< F(y, u{y), q, B 2 ) - F(x, u(x) ,p,B 2 + K 2 (Y, 5, rj)Id) , 

and, using the fact that Bi — B 2 < K 2 (Y,5,r])Id, we apply (H3b) to the left-hand side 
and (H2) to the right-hand side of (JSTj) . This yields 

XTi(B 2 - B 1 + K 2 (Y, 5, V )Id) < oji(\x - y\(l + \p\ + \q\) + w(\p\ A \q\)\p - q\)\\B\\ 

+ u 2 (K 2 (Y,5, V )) + C 1 

+ C 2 (\p\ 2 + \q\ 2 ) + C 3 \x - y\(\p\ 3 + \q\ 3 ) . (52) 

The estimates on the two sides of (|52p are done in the same way as in the proof of 
Theorem 12.11 The only difference is a term of the form CaO{rf)\Y\ a ~ 2 in the left-hand 
side of the estimate. 

Taking in account this additional term, we are lead to an analogous estimate to (|33l) 
with a right hand side of the form XK^ 1 + 0(r]) instead of \K~ l . We conclude in the 
same way by choosing first rj small enough. □ 

Step 2 : The case when G is independent of u but with a general dependence 
in p 

As for the treatment of the dependence in u, we are going to introduce a new variable. 
More precisely we introduce the function »:OxI->I defined by 

v(x,y) := u(x) - y . 

This new function is formally a solution of 

() V + F(x, v + y, D x v, D 2 xx v) = in O x E 



dy 2 
and 



-D,„vG I x, 1=0 in dO x 



DyV 



In fact, in order to justify this, one has just to be a little bit more precise about the 
definition of the boundary condition. We set, for x G dO, p x G M™ and p y G R 



G(x,p x ,p y ) 




if Py < 0, 
if Py > 0. 
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With this notation, the boundary condition for v becomes G(x, D x v, D y v) = and, be- 
cause in particular of the assumptions (G3)-(G4), it is rather easy to show that G, in 
addition to be homogeneous of degree 1 in (p x ,Py), satifies (G1)-(G2). 

On the other hand, the assumptions on the equation can be checked easily and there- 
fore the conclusion follows from Step 1. 
Step 3 : The general case 

In order to treat the dependence in u, as mentioned above, we use the method of the 
second step of the proof of Theorem 12.11 We are not going to give all the details since 
they are essentially the same. We just want to point out that in order to take care of 
the dependence in u and to be sure that the transformed boundary condition actually 
satisfies (G3)-(G4), one has first to introduce the function Gr defined by R > large 
enough and for x G dO, «6R and p G M n , by 

{G(x, u,p) if \u\ < R, 
G(x,—R,p) if u < — R, 
G(x,R,p) if u>R. 

Clearly, if R is large enough, u is still a solution of the Neumann problem with Gr and 
this transformation prevents difficulties with the behavior of G in u for \u\ large. 

And the proof of Theorem 12.21 is complete. □ 



Remark 2.2 We want to point out that, in (H2) ; the w term is needed only because we 
want to obtain local estimates : this is clear in the proof since, in ( fffffj) . this w term is 
used to take care of Q2 which comes from the localization term. Therefore, in the case 
of global estimates, the same result holds with w = 1. This remark can be used either 
on bounded domains or in unbounded domains where, if the equation and the boundary 
condition satisfy suitable uniform properties, L can be taken as small as we want (a u mild" 
localization) and the same effect occurs. 

3 The construction of the test-functions 

In this section, we provide the proof of Lemma 12.21 and 12.61 In particular we show how 
we construct the functions 0o and G5 which are used in the proof of Lemma 12.11 and 
12. 51 We will consider separately as in Section 1 the case of linear and nonlinear boundary 
conditions. 

3.1 The test-function for linear boundary conditions 

In this subsection we consider the case 

Ou 

— — h g(x) = on dO , 
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where 7: 80 — > M. n is a locally Lipschitz continuous vector field such that (7(2;), n(x)) > 
f > for any x G dO, and g: 90 — > M n is either a locally Lipschitz continuous or a 
locally Holder continuous scalar function. 

According to assumption (H3a), there exists a function A(-) G C^ c (dO, S n ) such 
that, for any x G 90, n(x) = A(x)j(x) and A (a;) > cold, for some constant cq > and 
such that ([9]) holds. Of course, this last property is the most important information in 
(H3a), the existence of such A without the connection with the ellipticity of the equation, 
being easy to show. 

We may assume without loss of generality, that ( r y(x),n(x)) = 1 for any x G dO, 

otherwise we change 7 in (7( jff (a;)) , 9 in M °)*n(x)) and A ( x ) m (l( x )i n ( x )) A ( x ); tnese 
transformations do not change the properties of 7 and g. 

Proof of Lemma 12. 21 As in the proof of the comparison result for this kind of problems 
(cf. Barles [3]), we are going to use regularizations of A and g. To do so, it is convenient 
to introduce the following lemma whose proof is classical and therefore left to the reader. 

Lemma 3.1 Assume that p G D(R n ), p>0, supp(p) C -8(0,1) and J Rn p(y) dy = 1. // 
/ G O 0,/3 (lR n ) for some < < 1, and f is bounded, then the function f : R n x [0, +00) — > 
R defined, by 

/>,£):=/ /(z)p(£^£)ldz for x G M™, e > 0, 

f(x,0) = f(x) forx£R n , 
is in O 0,/3 (lR n x [0, +00)). Moreover, the function f is C 2 in W 1 x (0, +00) with 

\DlJ{x,e)\, \DlJ{x,e)\, \D 2 J(x,e)\ < Ke^ 2 in K™ x (0, +00) 
for some constant K depending only on p, the L°° and the Holder norm of f . 

Step 1. The functions A and g and their regularizations. 

We first extend g and A to M n ; we still denote by g and A these extensions. We 
may assume that these extensions are respectively in in O ' /3 (R ?l ) and O ' 1 (R ?1 ). For some 
function p satisfying the properties of Lemma 13.11 (which is chosen and fixed from now 
on), we consider the functions A and g associated to A and g as in this lemma. Finally, 
we introduce, for some 5 > 0, the following quantity which is defined for £ G M. n by 

and we set 

A s (x, = A(x, A 5 (0) , g S (x, = g(x, A s (0) ■ 
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According to Lemma [3.11 these functions are C 2 as long as 5 > 0. 
We also observe that 

|AA*(£)|<1, |L>f e A a (£)| < KAJ 1 . 

Step 2. Construction of the functions ipo, ips and their main properties. 
For 5 > 0, we introduce the following function, for X, Y e R n and T > 

MX, Y, T) = (A 5 (X, Y)Y, Y) + K 1 5(2M - T) , (53) 

where K\ > is a constant to be chosen later and M is chosen so that 2M — T remains 
bounded. Moreover we set, for x and y in a suitable neighborhood of x Q 

if) S {x, y) := V>« ( — ^— , x - y, tZ(s) + %) 

We observe that, as 5 — > 0, ips an d converge locally uniformly respectively to ipo 
and ipo. Depending on the simplicity, we provide below result either on ips or ips, the 
translation from one to the other being straightforward. Most of the time we will use ips- 

In the sequel K > denotes a constant which may vary from line to line but which 
depends only on the data of the problem and is independent of the small parameter 5. 

Proposition 3.1 We have, for any X, Y e R n 

K^\Y\ 2 < MX,Y) < K\Y\ 2 } (54) 
K^\Y\ 2 < 4>s(X, Y) < K\Y\ 2 + K8 . (55) 

The proposition is straightforward consequence of the fact that A(x) > CqIoI for all 
i£l" Next we examine the regularity properties and the estimates on tpg and ips- 

Proposition 3.2 We have, for any X, Y G R n , T E R 

\D Y MX,Y,T)\<K\Y\ , \D X MX,Y,T)\<K\Y\ 2 , 
\D 2 YY MX,Y,T)\ < K , \D 2 XX MX,Y,T)\<K\Y\, 
\D 2 XY MX,Y,T)\<K\Y\, 

\D T $ S {X, Y,T)\<K5, D 2 TT ij s = D 2 TX ^ S = D 2 TY ^ 5 = . 

Moreover 

(DyMX,Y,T),Y) = 24> 5 + 0(\Y\ 3 ) + 0(8), 
(D 2 YY ^ S Y,Y) = 2^ s + 0(\Y\ 3 )+0(6), 

as Y — > and 5^0. 
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The proof of these estimates is tedious but straightforward: the main reason to provide 
Lemma [3.11 and to write A$ with a dependence in x and A$(£) was to have a simple way 
to check these computations. 
Proof of Proposition 13.21 We have 

D X MX,Y,T) = (D X A S Y,Y), 

D Y MX,Y,T) = (D Y AsY, Y) + 2A&Y , 

D 2 XX MX,Y,T) = (D XX A S Y,Y), 

D 2 YY MX,Y,T) = (D YY A S Y,Y) + 2A 5 + 2D Y A S Y. 

We premise some useful estimates on the first and second derivatives of Ag. By using 
Lemma [3.11 and the estimates on the first and second derivatives of A,5, we have 

D x A s (X,Y)=D As A = Y (l), 
D xx As(X,Y) = D 2 As A = AJ 1 O y (1), 

D Y A S (X, Y) = D As A^ = Y (1) , 

D YY As(X, Y) = (D As A)^ + Dl s A ^ = AfO Y (l) . 

Using the estimates on the first and second derivatives of properties of Ag, we obtain 
easily all the estimates on the derivatives of tps and the second part of Proposition 13. 21 □ 

We turn to the properties of ips with respect to the boundary condition. 

Proposition 3.3 // \x — y\ is small enough and K 1 is large enough, then we have, for 
some K > 

(D x xJ,s(x,y)Mx)) > -K\x-y\ 2 if x E dO , (56) 

(-D y M^y)n(y))<K\x- y \ 2 ifyedO. (57) 

Proof of Proposition 13.31 . We only check ( |56l) the other case being similar. 
By a direct computation, we have 

D x ^ s {x : y) = 2A 5 (x-y) + ((^^ + D Y A s )Y,Y)-K 1 5Dd(x) 
= 2A 5 (x-y)+0{\x-y\ 2 ) + K 1 5n(x) 

But we recall that 

\\A s -A\\<K(\x-y\ + 5) 
and since x G dO, by the regularity of the boundary, we have 

(n(x), (x — y)) = d(y) + 0(\x — y\ 2 ) as \x — y\ — > . 
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Thus if Ki is large enough we have 

(D x iPs(x,y),l{x)) > (2A(x-y), 1 )-K\x-y\(\x-y\+8)-K\x-y\ 2 + K 1 8, 

> (2(x-y),n(x))-K5-K\x-y\ 2 + K l 6, 

> -K\x-y\ 2 . 

□ 

Step 3. Construction of the functions Xo an d Xs- 
In the same way as above, we set for 5 > 

Xs(X, Y, T, Z) = g s (X, Y)Z + K 2 S (2M - T) , (58) 

where K-i > is a constant to be chosen later and M > is chosen as above. We also set 

Xs(x, y) := Xs ( ^-tt-, x-y, d(x) + d(y), d(x) - d(y) 



One can easily check that xs and Xs converges locally uniformly respectively to xo an d 
Xo as 5 — > 0. 

As in the previous step, we first consider the regularity properties of xs- 



Proposition 3.4 For every 5 > we have 



DxXs(X, Y, T, Z)| < ATA?" 1 ^ , \D Y xs(X,Y,T,Z)\ < AAf' : 



\D 2 xs(X,Y,T,Z)\<KA p s - 2 \Z\ . 

Proof of Proposition 13.41 Again the tedious computations are simplified by the way 
we write down g$. We first observe that 

D x g 5 = D x ~g = kl l O Y {\), 
D 2 xx ~g 5 = D xx g = A^ 2 Y (1) , 

DrSs = D As g = Af-'Oyfl) , 

^ = Au(^) 2 +c Ai ^=Aro r( i). 

Then, by the definition of xs, we deduce 

D x xs(X,Y,T,Z) = D x g 5 Z<KA P s - l \Z\, 

D Y xs(X,Y,T,Z) = D Y ~g & Z<KhP & - x \Z\, 

D 2 XX UX,Y,T,Z) = D 2 xx ~g 5 Z < KA^ 2 \Z\ , 

DyyXs(X, Y, T, Z) = D 2 YY ~g 5 Z <KA^ 2 \Z\. 



□ 



Then we consider the boundary condition. 
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Proposition 3.5 If \x — y\ is small enough and K 2 is large enough, then we have, for 
some K > 

(D x X6(x,y), 1 (x))+g(x)> -K\x-yf + K5? if x e dO , (59) 
(-D yX s(x,y), 1 (y))+g(y)<K\x-yf-K5P ifyedO. (60) 

Proof of Proposition [375] . Again we only check the first property (1591) . We first notice 
that, by the definition and properties of gs, 

IH^ 1 , x-y)- g(x)\ <K(\x-yf + 
On the other hand, we have 

D x xs(x,y) = g 5 Dd(x) + {^D x gs + D Y gs){d{x) - d(y)) - K 2 5 l3 Dd(x). 

Therefore since \x — y\ is small, we obtain for K 2 large enough 

x ~\~ y 

(D x xs(x,y),j(x)) + g(x) > [g(x) - 9s(—^— , x - y)\ 

+ ((-D x g s + D Y g s ),j(x))(d(x)-d(y)) 
+ K 2 5? 

> -K\x - yf + KS P , 

Thus we have shown (1591) . □ 
Now we can prove the following result. 

Proposition 3.6 If \x — y\ << 1, then we have 

(D x e s (x,y), 1 (x))+g(x)>0 if x e 80 , (61) 
(-D y e s (x,y)^(y)) + g(y)<0 ifyedO. (62) 

Proof of Proposition [3761 Again we only check the first property (I6TI) . We first observe 
that, because of the properties of ips and since d(x) = 0, we have 

(D x (e-^ +d ^(Mx,y))%l(x)) = Ke~ Rd ^Mx,y) 

+ | e -^%)(^ (X; y))^)( Dx Mx, V),l(x)) 

> K e -Kd[y) | x _ y] a _ Ke -Kd(y) \ x _ y |« 

Thus, if K is sufficiently large, we obtain 



(A, 



-KW*Mv))ty s ( Xi y))S] ;7 ( x )) > R\x-y\ 
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for some constant K. Similarly one can shows that 

(D x [e-^)^ ))^^,^))] , 7 (a0> > K\x - x \ 2 . 

Then if \x — y\ « 1 and C > is large enough, we have 

(D x Q s (x,y), 1 (x))+g(x) = {CD x [e- R ^ +d ^\U^y))H^)) 

+ (LD x [e- kd ^Mx,x )] 7 j(x)) 

+ (D x xs(x,y) } -f) + g(x) 

> CKe- R{d{x)+d{y)) \x-y\ a -K\x-yf + K5^>0. L 

The proof of Lemma 12.21 is obtained by combining Propositions 13.21 13.41 and 13.61 
We conclude this section with the following lemma. 

Lemma 3.2 If P = A^ 1 (x)D x Qs(x,y) , we have 

p = Y + oy(l) + 0(8) as \Y\ ^0 and 5^0. 
Proof. We just give a sketch of proof. We first observe that 

\\M^y, x-y)- A(x) 1 1 <K(\x-y\+5). 

Therefore, by direct computations, we obtain 

D x Q s (x,y) = Ce- R ^ x)+d ^ (v^)^ [AY + 0(\Y\ 2 ) + O 
= Ce- R{d{x)+d{y)) te)^ 21 [AY + 0(\Y\ 2 ) + +0 
Thus we get 

A-\x)D x Q s {x,y) Y + Q(\Y\ 2 ) + Q(S) 

\A-i(x)D x Q s (x,y)\ \Y + 0(\Y\ 2 ) + 0(5)\ 

= Y + o y (l) +0(5) as \Y\ -> 0,5^ 0, 

which gives the desired result. 

3.2 The test-function for nonlinear boundary conditions 

We recall that we have to build this test-function in the case when the function G 
independent of u and homogeneous of degree 1 with respect to p. 
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We first extend the function G(x,p) to M. n x R n and we may assume that all the 
properties of G are still satisfied in V x R" where V is a neighborhood of dO. The 
properties of G imply that, for every x G V, p G M n there exists a unique solution 
t = (p(x,p) of the equation 

G(x,p + tn(x)) = 0. (63) 

One can verify that <p is still homogeneous of degree 1 and satisfies (G2). 

It is not restrictive to reduce to the case when the boundary is flat and more precisely 
= {x n > 0} and dO = {x n = 0}. 

Proof of Lemma 12.21 

Step 1. The function ip and its regularization. 

In order to regularize the function ip, we first extend it to MJ 1 and we still denote by ip 
this extension. We may assume that this extension satisfies (G2). 

We introduce, for S, rj > 0, the following quantity which is defined for £ G M. n by 

r, 



r* (0 = v 

and we set, for x, £ G 

<ps(x,£) = tp(x,£,T s (£)) , 

where <p is defined as in Lemma 13.11 and £ = 4r . 

We first observe that the following estimates, which are used extensively in the sequel, 
hold 

Krj Krj 

\ D ^\ < (1+1^2)1/2 » l^a r «l < (l + |£|2)l/2 (1 + | e j2 ) l/2 • 

Step 2. Construction of the functions ipo, ipg, and their main properties. 

For 5 > 0, we introduce the following function, for X, Y G M n with as above Y = 
(Yi,...,y n )andX = (X 1 ,...,X n ) 

MX, Y) = \Y\ 2 - 2^ 5 (X, Y)\Y\Y n + 2A 1 Y^ + K^R - X n ) , (64) 

with Ai,K\ > constants to be chosen later. The constant R > has to be chosen in 
order that the term R — X n remains positive; this does not create any problem since we 
argue locally. Moreover we set, for x and y in a suitable neighborhood of Xq 

ip s (x, y) := tps ( — y~ ' x ~ V 

We observe that, as 5 — ► 0, ips and ^<5 converge locally uniformly respectively to an d 
ipQ. As in previous subsection we provide below result either on ip§ or ip§, the translation 
from one to the other being straightforward. Most of the time we will use ipg. 

In the sequel K > will denote a nonnegative constant which may vary from line to 
line but which depends only on the data of the problem and is independent of the small 
parameters 5 and r\ 
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Proposition 3.7 If Ai > is large enough, we have, for any X, Y G M n 

K^\Y\ 2 < $ (X,Y) < K\Y\ 2 } (65) 
K^\Y\ 2 < 4> S (X,Y) < K\Y\ 2 + K5. (66) 

We skip the proof of this proposition which is straightforward : it is based only on 
Cauchy-Schwarz's inequality to control the (ps or the (po term and on the fact that (ps and 
(fo are bounded. 

Next we examine the regularity properties and the estimates on ips and ips- 

Proposition 3.8 If the constant r] > is chosen small enough and Ax , K\ large enough, 
then, for any 5 small enough and for all X, Y G IR n , we have 

\DxMX,Y)\ < K\Y\ 2 + K5 , \D Y ips(X,Y)\ < K\Y\ , 

(D Y MX, Y),Y) = 2MX, Y) + 0(<n)0(\Y\ 2 ) + 0(6) , 
(D 2 YY MX, Y)Y, Y) = 2MX, Y) + 0( V )0(\Y\ 2 ) + 0(6) , 
as Y — > 0, i] — > and 6^0, 

\D 2 XX MX,Y)\ < j\Y\(6 2 + |Yf) 1/2 , \D 2 XY MX,Y)\ < K\Y\ , \D 2 YY iPs(X,Y)\ < K. 

Proof. The proof of these estimates is tedious but straightforward. Lemma 13.11 and 
the way we write (p with a dependence in x, £ and Ts(£) is a simple way to check these 
computations. 

By direct computations, we have 

Dx$ s (X,Y) = -2D xl p 5 \Y\Y n -K5e n , 
D Y iPs(X,Y) = 2Y - D Y ip s \Y\Y n 

Y . , 

~ 2t fiSTyi Y n - 2(ps\Y\e n + 2A l Y n e n , 

D 2 XX MX,Y) = -2D xxl f S \Y\Y n , 
D 2 YY MX,Y) = 2Id-2D YY <ps\Y\Y n 

Y . . 

-AD Y (p s —Y n - AD Y ips\Y\e n , 

Y Y 

-4(^ 5 — e n - A^sD Y — Y n + 2A{Y n , 

where e n is the n-th vector the canonical basis of IR n . In this case, we also have n(x) = —e n , 
for all x G dO. 
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We estimate the first and second derivatives of (pg. By using Lemma 13.11 and the 
estimates on the first and second derivatives of T$, we have 

\D x s (X,Y)\ = \D x <f(X,^- T S {Y))\<K, 

\D 2 xx 6 (X,Y)\ = \D XX 0(X,^V S (Y))\<^-, 

\D Y s (X,Y)\ = \Ds<pD Y — + D c 0D Y r s \ < — 



\D YY <p s (X,Y) = \DI&(dy^ +D i tpD YY ^ + 2Dl c {pD Y ^®D Y T s 



+D 2 (C (D Y T s y + D c <pD YY T s \ 
K A 

By combining the above estimates we obtain 

\D X MX,Y)\ < K\Y\ 2 + 0(5), \D Y MX,Y)\<K\Y\, 

\D 2 XX MX,Y)\ <r / - 1 iv|F|(|F| 2 + 5 2 ) 1 / 2 , \D Y tjjs(X, Y)\ < K\Y\ . 

Next we estimate (Dyif)g(X,Y),Y) and {D YY ifjgY, Y ) . Tedious but straightforward com- 
putations show that 

I (D Y T 5 , Y}\ < Krj, | (D YY T 5 Y, Y)\ < K V . (67) 

Moreover 

(Dy4s(X,Y),Y) = 2\Y\ 2 -4<p 5 \Y\Y n + 2A 1 Y 2 

-D C if S ((D Y T 5 ,Y))\Y\Y n 
= 2^ 5 + 0(5) + 0(r ] )\Y\ 2 ; 



(D YY i> s Y,Y) = 2\Y\ 2 -A0 s \Y\Y n + 2A 1 Y 2 

-2D ( if S (D 2 T 5 Y,Y)\Y\Y n -AD^ s ((DT s ,Y))\Y\Y n 
-Dl^ 5 ((DT s ,Y)) 2 \Y\Y n 

= 2$ s + 0(5) + 0(r))\Y\ 2 . 

These properties complete the proof of Proposition 13.81 □ 
We turn to the properties of ips with respect to the boundary condition. 
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Proposition 3.9 If \x — y\ is small enough and Ki is large enough, then we have, for 
some K > 



G{x, D x ip s (x, y)) > —K\x - y\ 2 if x n = , 
G(y, -D v if> s (x,y)) < K\x - y\ 2 ify n = 0. 



(68) 
(69) 



Proof. We only check (1681) the other case being similar. 

If x G dO then d(x) = x n = 0; moreover since y G O, d(y) — y n > 0. Thus we have 

x ~\~ y 

D x ip s {x,y) = 2{x-y)-2{p s {——,x-y)\x-y\e n 

+ (d x 0s(^-^-, x-y) + D Y (fs{^-^-, x - y)J \x - y\y n 
- 2ip & {-— — ,x - y)- -y n -2A x y n e n -K 1 6e n . 

2 \x-y\ 

By taking in account that n(x) = — e n we have 

D x i/j s (x,y) =p + q + r, 



where 

p = 2(x — y) + 2<p(x, x — y)n{x) , 

Q = 2 [0s{^Y^-,x-y)\x-y\-cp(x,x-y)]n(x) 

+ (D x ^s(^—^-,x- y) + D Y <fs{^-^ y))\x - y\y n 
~ ( x + y ^x-y 
2 \x-y\ 
r = (2A 1 y n + K 1 5)n(x). 

We first notice that, taking in account the definition of fig, since <p is homogeneous of 
degree 1 with respect to p and satisfies (G2) we have 



~ t x+ y ^ i x ~ y m i 

<P6(—z—> x -y) -<Pfa] r ) x-y 

2 \x- y\ 



< K\x-y\(\x-y\ + F s ) 

< K[\x - y\ 2 + 5 + (y n - x n )\ 

< K\x-y\ 2 + K5 + Ky n . 



Moreover, we have 



x - y) l_ v y \ yn\ < K \ x - y\ 2 - 
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By using the fact that G(x,p) = and combining the above estimates with the properties 
of G, we get 

G(x,Djj> s ) > G(x,p + r) - K\q\ 

> \{2A 1 y n + K 1 5)-K\x-y\ 2 -K5-Ky n . 

Thus if Ai and K\ are large enough we get 

G(x, D x ^ 5 ) > -K\x - y\ 2 + Kb . 

□ 

Step 4. Construction of the function 0o and its main properties. 
For C, A2 > large enough, we introduce the following function 

o (x, y) = \ x - x \ 4 - 2C\x - x \ 3 (d(x) - d(x )) + A 2 (d(x) - d(x )) 4 . (70) 

In order to check that the function 0$ satisfy the right boundary conditions we premise 
the following Lemma whose proof can be found in [5]. To formulate it, we use the following 
notation : for p G M n and x G dO, T(p) := p — (p,n(x))n(x). T(p) represents the 
projection of p on the tangent hyperplane to 80 at x. 

Lemma 3.3 Assume that (Gl) and (G2) hold and that, for some x G dO, and p G W 1 , 

we have G(x,p) < (resp. G(x,p) > 0), then there exists a constant K (depending on v 
and K in (Gl)-(G2) y ) such that, if (p,n(x)) < —K\T(p)\, then 

G(x,p + p) < , 

(resp. if (p,n(x)) > K\T{p)\, then 

G(x,p + p) > ). 
The connection with 0o is given by the following result. 
Lemma 3.4 For C large enough and for all x G dO, <po satisfies 

^(x) > K\T(D<f> )\ , 

where K > is the constant appearing in Lemma \3. & 

Proof. We first compute the normal derivative of <fio and we use the usual property 
linking the distance function and n; this yields 

(D x (f) (x),n(x)) = (4\x — xo\ 2 (x — xo),n(x)) + 2C\x — xo\ 3 

= 4\x - x \ 2 (d(x ) + 0(\x - Xq\ 2 )) + 2C\x - x | 3 

> 2C\x - x \ 3 + 0(\x - x \ 4 ) . 
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On an other hand, we clearly have \T(D(p )\ < 4| X — Xq | 3 . Thus, if K is the constant 
given in Lemma [3.31 by choosing C large enough and x close to Xq, we have 



2C\x — Xq 


3 + 0(\x — x \ 


4 ) 


4| 


X — Xq\ 


3 



\T(D ( j )0 )\- L (D x( j )0 (x),n(x)) > — ! £ ^ ^>K 



□ 

Now we can prove the following result. 
Proposition 3.10 If \x — y\ << 1, then we have 

G(x,D x e 5 (x,y))>0 ifx n = 0, (71) 
G(y,-D v e s (x,y))<0 ify n = 0. (72) 

Proof of Proposition [37101 Again we only check the first property (1711) . First of all we 
observe that, because of the assumption (Gl) and the property (1681) of ips we have 

G(x,D x (e- R{y "\Mx,y)) f ) = G{x,Ke- k ^$ s {x,y))?n(x) 

+ ^e-^(iP s (x,y))%- l D x ^ s (x,y)) 

> fike- RiVn) \x - y\ a - Ke- R(yn) \x - y\ a . 

Thus, if K is sufficiently large, we obtain 

G(x, D x (e- Rix " +y "\Mx, y)) § )> K\x - y\ a , 

for some constant K. 

Now by combining Lemma I3.3I and I3.4I we get 

G{x,D x Q s {x,y)) = G^CD^e-^^iM^y^ + LD^o) 
> K\x-y\ a . 

And the result is proved. □ 
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